arXivil 504.07055vl [cond-mat.mes-hall] 27 Apr 2015 


High-Q Optical Cavities in Hyperuniform Disordered Materials 


Timothy Amoah and Marian FlorescrQ 
Advanced Technology Institute and Department of Physics, 

University of Surrey, Guildford, GU2 7XH, United Kingdom 
(Dated: April 28, 2015) 

We introduce the first designs for high-Q photonic cavities in slab architectures in hyperuniform 
disordered solids displaying isotropic band gaps. Despite their disordered character, hyperuniform 
disordered structures have the ability to tightly confine the TE-polarised radiation in slab config¬ 
urations that are readily fabricable. The architectures are based on carefully designed local modi¬ 
fications of otherwise unperturbed hyperuniform dielectric structures. We identify a wide range of 
confined cavity modes, which can be classified according to their approximate symmetry (monopole, 
dipole, quadrupole, etc.) of the confined electromagnetic wave pattern. We demonstrate that qual¬ 
ity factors (Q) Q > 10 9 can be achieved for purely 2D structures, and that for three-dimensional 
finite-height photonic slabs, quality factors Q > 20, 000 can be maintained. 

PACS numbers: 41.20.Jb, 42.70.Qs, 78.66.Vs, 85.60.Jb 


A special class of disordered photonic heterostructures 
has recently been shown to display large isotropic band 
gaps comparable in width to band gaps found in pho¬ 
tonic crystals m- The large band gaps found in these 
structures are facilitated by the hyperuniform geometri¬ 
cal properties of the underlying point-pattern template 
upon which the structures are built. The statistical 
isotropy of the photonic properties of these materials is 
highly relevant for a series of novel photonic functional¬ 
ities including arbitrary angle emission/absorption and 
free-form wave-guiding 0 , a. 

A point pattern in real space is hyperuniform if for 
large R the number variance cr(R) 2 within a spherical 
sampling window of radius R (in d dimensions), grows 
more slowly than the window volume, i.e., more slowly 
than R d . In Fourier space, hyperuniformity means that 
the structure factor S( k) approaches zero as |k| -A 0. 
016 ]. The lack of periodicity in these hyperuniform dis¬ 
ordered solids (HUDS) demonstrates that Bragg scatter¬ 
ing is not a prerequisite for photonic band gaps and that 
interactions between local resonances and multiple scat¬ 
tering are sufficient, provided that the disorder is con¬ 
strained to be hyperuniform [I]. 

The concept of optical cavities in hyperuniform dis¬ 
ordered (HUD) photonic materials was recently intro¬ 
duced in Ref. [7]. The structure analysed was obtained 
by placing dielectric rods at each point of a hyperuniform 
point pattern. The radius of one selected rod was varied 
to achieve localisation of the TM-polarised electromag¬ 
netic field at that point. The study was based purely 
on 2D structures and vertical confinement, the primary 
loss pathway in real slab structures, was not discussed. 
The question of how to achieve index-guiding and the es¬ 
sential vertical confinement in disordered photonic slab 
structures, a prerequisite of realising cavities with high- 
quality factors in fabricable structures, was seen as a po- 
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tential roadblock in the HUD photonic materials field. 

In this paper we discuss analogous network structures 
for localising TE-polarised radiation and analyse the ver¬ 
tical confinement issue. The structures analysed are com¬ 
posed of polygonal cells of dielectric walls. The protocol 
for generating these structures described in Refs. mm 
consists of Delaunay triangulating over a hyperuniform 
point pattern and connecting the center of mass of the 
Delaunay triangles to form cells that contain one original 
point each. The walls of the resulting network are then 
given a finite thickness. Trivalency is preserved through¬ 
out the structure. The structure presents well defined 
short-range order and long-range statistical isotropy. The 
resulting structure can be described as a hyperuniform 
disordered honeycomb (HDH). Recently, these designs 
have been fabricated on the microwave scale and suc¬ 
cessfully tested [9, TO] . 


FIG. 1: (Color online) Magnetic field distributions of two 
Anderson-like localised modes, (left) below the lower and 
(right) above the upper photonic band gap edge respectively. 

We begin our analysis with the study of purely 2D 
structures. A length scale a = L/y/N is defined, such 
that an 7V-point hyperuniform pattern in a square box 
of side length L has density of 1/a 2 . The structure pa¬ 
rameters are set to e = 11.56 for the dielectric constant 
and the wall width is set to w = 0.2 a. For our study, 
we employ a point pattern containing N = 500, with a 
hyper uniformity order parameter x = 0.5 m- We cal¬ 
culate the photonic band structure using the eigenmode 
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expansion software “MPB” for a y/Na x y/Na sample 
and identify the photonic band edges. We then run a 
sweep for varying wall thickness (w) and find the largest 
band gap with Auj/ujc = 30.7% for w = 0.23a, where 
Auj is the band gap width and ujc is the band gap cen¬ 
ter frequency. Similar to other conventional disordered 
structures, unperturbed HUD structures display defect 
modes. These Anderson-like localised modes, shown in 
Figure [lj occur naturally at the photonic band gap edges 
and extend over over network domains comprising 5-10 
cells. Sparsely occurring are accidental localised modes 
extending on 1-3 network cells which are promoted into 
the photonic band gap topologically, but we are not ex¬ 
ploring these here [12]. 

In an otherwise unperturbed HDH structure, it is pos¬ 
sible to create an intentional localized state of the elec¬ 
tromagnetic field by reducing or enhancing the dielectric 
constant at a certain point in the structure. For a trian¬ 
gular lattice of holes it is common practice to fill a single 
hole to make a cavity which is often labelled a HI cav¬ 
ity. In analogy to this we fill a single cell and also label 
it HI. Due to the presence of the defect, four localized 
cavity modes are created within the photonic band gap 
at specific frequencies. The mode profiles are shown in 
Figure [2] Two of the modes are dipole-like and the other 
two modes are quadrupole-like. We denote the lower fre¬ 
quency dipole mode Di and the higher frequency dipole 
mode D 2 . The cavity has an approximately hexagonal 
shape, as such we can describe the modes according to 
the approximate symmetries associated with a hexago¬ 
nal unit cell in a triangular lattice. For the Di mode the 
nodal axis of the H z component of the electromagnetic 
field lies approximately along the faux K direction and 
the mode is propagation confined along the faux M direc¬ 
tion. For the D 2 mode the nodal axis of the H z compo¬ 
nent of the electromagnetic field lies approximately along 
the faux M direction and the mode is propagation con¬ 
fined along the faux K direction. The magnetic field dis¬ 
tribution of the cavity modes shown in Figure [2] suggests 
very good in-plane confinement and calculations of the 
quality factors (using the finite difference time domain 
software “MEEP” [13]) confirms that the quality factor 
for all modes higher than 10 9 . 


Mode 

uj/ujq in-plane 

uj/ujq slab 

Q-factor 

Di 

0.22482 

0.26722 

7976 

d 2 

0.23775 

0.28003 

8672 

Q 

0.24546 

0.28811 

2561 

H 

0.25133 

0.29385 

4845 

O 

0.25916 

0.29973 

3230 


TABLE I: The mode frequencies in the 2D and 3D case and 
Q factor for the modified HI cavity with a 0.2a inner hole 

To study the frequency behaviour of the cavity modes 
as a function of dielectric filling of the defect cell, we are 
placing a circular air hole at the center of the unmodi¬ 
fied HI cavity, vary its radius from r = 0 to r = 0.535a 
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FIG. 2: (Color online) Magnetic field distribution for the 
four defect modes in the unmodified HI cavity in a structure 
of wall thickness w — 0.2 a. The modes are labelled Di, D 2 , 
Qi, Q 2 (shown from left to right in the figure). 
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FIG. 3: (Color online) Frequencies of the cavity modes in 
the band gap obtained for different radii of a hole placed at 
the center of the cavity. 


and monitor the evolution of the cavity mode frequen¬ 
cies. The general trend is that the frequency of the cav¬ 
ity modes increases with the radius of the hole. In other 
words, the mode frequency is lowered by the inclusion 
of more dielectric, which demonstrates that the cavity 
originates from modes above the upper band edge of the 
uperturbed structure. The same behaviour is observed in 
the TM case when the radius of a defect rod is increased 
0 . 

We now turn our investigation towards the photonic 
slab where we no longer can neglect the vertical ex¬ 
tent of the structure. This is the same as the concept 
of photonic-crystal slabs [14] [T5] . which has proven to 
be a promising platform for integrated photonic micro¬ 
circuitry. In photonic slabs, while there exists no true 
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FIG. 4: (Color online) Magnetic field distribution for the five cavity modes in the enlarged HI cavity with center hole r = 0.2a 
in a structure of wall thickness w m. 0.4a . The modes are labelled Di, D 2 , Q, H, 0 left to right. 



FIG. 5: (Color online) Frequencies of the cavity modes in the 
band gap obtained for different radii of a hole placed at the 
center of the cavity. Modes are labelled according to their 
number of pole lobes. Different colours are used for modes 
with the same number of lobes but which have distinct field 
patterns. Cross over intermediate modes are coloured black. 
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FIG. 6: (Color online) Intensity distribution (yellow) of a 
confined cavity mode in a hyperuniform disordered honey¬ 
comb (HDH) photonic slab (blue). 


band gap, the existence of “pseudo-band-gaps” enables 
low-loss waveguiding and high Q cavities [17]. The 
“pseudo-band-gap” is associated with an effective band 
gap in the in-plane projected band structure for which 
modes below the light cone cannot couple to the contin¬ 


uum of states outside the slab. This can also be thought 
of in terms of “index confinement” due to the contrast in 
effective dielectric inside the slab and in the surrounding 
air region. The projected band structure method only 
works well when only a single in-plane unit cell is con¬ 
sidered. For supercell calculations, the spectral regions 
with no effective band gap will be translated onto the 
region of the effective band gap, through the “supercell 
folding” effect. With no straightforward method to cir¬ 
cumvent this we will resort to making the assumptions 
that a pseudo-band-gap does exist for the photonic slab 
and that if a cavity mode of specific form lies within the 
band gap when calculated in 2D it will also lie within the 
pseudo-band-gap for the photonic slab. Consequently, 
the frequencies in Figure [3] are the in-plane frequencies, 
not the slab frequencies. 

For photonic slabs it is not generally true, that the 
best optical confinement is achieved for the largest in¬ 
plane photonic band gap. For a large enough sample the 
main losses occur in the vertical direction and to prevent 
them one needs to rely on index confinement. As such, 
the thicker walls translate into a higher effective index 
for the slab, and hence, better vertical confinement. For 
our structures, we use a structure height of h = 0.7a, and 
doubled the wall width to w = 0.4a to improve confine¬ 
ment. For an unmodified HI cavity the expected quality 
factor is typically rather moderate even for a conventional 
photonic crystal, Q < 500 [IB]. We calculate Q = 210 for 
the Di and Q = 190 for the D 2 mode in the HDH. This 
values can be considerably increased by placing an inner 
hole at the centrer of the cavity (see Table |I|. 

A popular method of of achieving optimal cavity de¬ 
signs is to reduce adjacent holes slightly in size and shift 
them outwards along the lattice directions. In the HDH 
case there are no lattice directions, and we instead shift 
the shrunk cells along the vector given by the center of 
mass of the cavity to the center of mass of the neigh¬ 
bouring cell. The cells are shrunk to 52% of the size of a 
cell with infinitesimal walls and are shifted by 8% of the 
length of the vector outwards. For a modified cell cav¬ 
ity we find 5 cavity modes laying in the band gap. The 
lowest frequency mode is again the Di mode followed by 
the D 2 mode. An isocontour plot of the electric field in¬ 
tensity distribution of the Di mode is shown in Figure [6] 
displaying in-plane PBG-mediated confinement and ver- 
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tical refractive index confinement. As expected, the two 
dipole modes have now lower frequencies due to the ad¬ 
dition of dielectric material. The next higher frequency 
modes are a quadrupole-like mode (Q) and a hexapole- 
like mode (H). Comparing the two mode patterns, we 
note that (due to disorder) the symmetry separation of 
the two modes is not entirely complete. Lastly, there is 
the highest frequency mode which is rather difficult to de¬ 
fine. The H z field is symmetric with respect to the faux 
K direction; therefore the mode is an odd mode. We can 
identify 4 nodes lying “within” the cavity, so it can either 
be considered a second quadrupole mode or an octapole 
mode. We assert that since it is higher in frequency than 
the hexapole mode we label it as an octapole-like mode 

(O). 

Figure [7] shows the spatially Fourier transformed in¬ 
plane fields of the Di and D 2 cavity modes for both the 
unmodified and modified (extra central hole) cavity. The 
respective light circles are indicated by the white dots. 
Only the Fourier components with k < ujc can radiate 
into the far-held and hence the more Fourier components 
that are within the light circle the greater the radiative 
loses in the vertical direction. Clearly, the modihed- 
cavity design in Figure [7] displays a cosniderably smaller 
number of radiative components. 


of mass shift of 8%, we find a very high quality factor of 
Q = 20,148, which is significantly larger than Q = 3230 
obtained at 52% neighbour cell size. Figure [8] shows the 
Fourier spectrum in each respective case. As before, the 
Fourier components in the light circle are reduced for 
the higher quality modification, however the difference 
in the reduction is less even if a large increase in qual¬ 
ity factor is achieved. Most significant is the change in 
frequency from 0.29973 to 0.28702, which displace the 
mode to spectral domain where the radiative losses are 
minimised. 




FIG. 7: (Color online) Spectrum of the Fourier components 
of the electric field distribution for the Di mode before (top 
left) and after modification (top right) and D 2 mode before 
(bottom left) and after modification (bottom right) . Loga¬ 
rithmic colour scale. 


Finally, we consider optimal cavity designs for the O- 
mode. For a neighbour cell size of 44% at the same center 


FIG. 8: (Color online) Spectrum of the Fourier components 
of the electric field distribution for a modified cavity with 
neighbour cell size of 52% (left) and 44% (right), both shifted 
by 8%. Logarithmic colour scale. 


In summary, we have introduced novel architectures 
for the design of optical cavities in a hyperuniform dis¬ 
ordered material. We have demonstrated that Hl-type 
cavity defects can support localized modes with a variety 
of symmetries and multiple frequencies. The ability to lo¬ 
calize modes of different symmetry and frequency in the 
same physical cavity and to guide light through modes 
with different localization properties can have a great 
impact on all-optical switching m, implementations of 
linear-optical quantum information processors [20 , and 
single-photon sources [21] 322] • 

In principle, one may have the presumption that disor¬ 
der would facilitate significant out of plane scattering as 
compared to periodic structures. Here we have demon¬ 
strated that for cavity built on hyperuniform platforms, 
by adequately adjusting the structure parameters and 
cavity design it is possible to achieve very tight optical 
confinement. 

Our successfully demonstration of high-Q cavity for 
transverse electric (TE) polarized radiation is encourag¬ 
ing future investigation of TE wave-guides in disordered 
photonic slabs. These would be promising candidates for 
achieving highly flexible and robust platforms for inte¬ 
grated optical micro-circuitry. 
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